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Resonant Electric Dipole-Dipole Interactions between Cold Rydberg Atoms in a 

Magnetic Field 
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Laser cooled *^Rb atoms were optically excited to 46d5/2 Rydberg states. A microwave pulse 
transferred a fraction of the atoms to the 47p3/2 Rydberg state. The resonant electric dipole-dipole 
interactions between atoms in these two states were probed using the linewidth of the two-photon 
microwave transition 46(^5/2 — 47^5/2- The presence of a weak magnetic field ~ 0.5 G reduced the 
observed line broadening, indicating that the interaction is suppressed by the field. The field removes 
some of the energy degeneracies responsible for the resonant interaction, and this is the basis for a 
quantitative model of the resulting suppression. A technique for the calibration of magnetic field 
strengths using the 34si/2 — 34pi/2 one-photon transition is also presented. 

PACS numbers: 32.80.Rm, 34.20.Cf, 32.80.Pj 



I. INTRODUCTION 

The interactions between neighboring gas phase atoms 
can be enhanced by exciting them to Rydberg states. 
To obtain translationally cold interacting Rydberg gases 
(« 100 ^iK), laser cooled atoms can be excited [l|, 12 ■ 
The behavior of these cold Rydberg gases is often quali- 
tatively different from that of hotter samples (w 300 K) . 
For example, a cold cloud of Rydberg atoms can sponta- 
neously evolve into a plasma under conditions in which 
a hot cloud could not 3]. Recent theoretical calcula- 
tions suggest that cold interacting Rydberg atoms on uni- 
formly spaced lattices will show properties analogous to 
that of a plasmon wave propagating along metal nanopar- 
ticles 3. 

Translationally cold Rydberg atoms have recently re- 
ceived attention in the context of neutral atom quantum 
information processing. Theoretical work has suggested 
that Rydberg atom interactions may be used to imple- 
ment gates between atoms holding qubits [3, for encod- 
ing qubits in multiatom clouds laK and to implement "on 
demand" single-photon sources [3 • Recent experimental 
work has demonstrated a suppression of the average Ryd- 
berg density due to interactions - so called "local block- 
ade" H, 13 J and rapid pro gres s is expected in this area 
over the next several years |lOj . 

Transitions between Rydberg states can be driven us- 
ing microwave sources, which offer stability, linewidth, 
and modulation capabilities unrivaled by lasers. To elu- 
cidate the role of multibody interactions in cold Rydberg 
gases, Mourachko et al. 11] used microwave radiation to 
transfer a fraction of the optically excited Rydberg atoms 
into a different angular momentum state. Recently, Li et 
al. [I3I have used a similar microwave redistribution pulse 
to establish the role of the electric dipole-dipole interac- 
tion in plasma formation from cold Rydberg gases. In 
addition, the linewidths of microwave driven transitions 
between Rydberg states can be used as a sensitive probe 
of interactions between Rydberg atoms [13 ■ 

Recently, we have probed the strong resonant electric 



dipole-dipole interactions between ^'^Rb atoms in A5d^/2 
and 46p3/2 states using the linewidth of a two-photon 
microwave transition from 45(^5/2 to 46(i5/2 [13 ■ The 
results were consistent with theoretical expectations, and 
this approach promises to be a powerful technique for 
studying cold Rydberg atom interactions. 

How will weak magnetic fields (w 0.5 G) influence the 
interactions between Rydberg atoms? Some of the en- 
ergy degeneracies which are responsible for the resonant 
electric dipole-dipole interaction observed in Ref. [13 can 
be spoiled by application of a magnetic field. It may be 
desirable to excite neutral atoms to Rydberg states in a 
loffe-Pritchard trap without switching off the trapping 
fields. Rydberg atoms themselves could also be trapped 
using inhomogeneous magnetic fields |l4| . 

In the present work we use microwave spectroscopy of 
Rydberg atoms to calibrate DC magnetic field strengths. 
Microwave fields are then used to manipulate Ryd- 
berg populations, and probe interactions, demonstrating 
that the dipole-dipole interaction between cold Rydberg 
atoms is partially suppressed by magnetic fields. A quan- 
titative model is presented which explains these results. 



II. EXPERIMENTAL 

The apparatus is similar to that described in Ref. [13 ■ 
Cold *^Rb atoms in a vapor-cell magneto-optical trap 
(MOT) [131 are excited to high-n Rydberg states using 
a pulse of laser light. Microwave pulses are then used 
to transfer Rydberg population and probe interactions. 
The final state populations are analyzed by selective field 
ionization (SFI) 16]. The significant differences between 
the present apparatus and that in Ref. 13] are the laser 
excitation source and the control of the magnetic field 
during the experiment. 

Instead of a pulsed dye laser system, a frequency- 
doubled cw Tiisapphire system is used to excite Ryd- 
berg states. A commercial Ti:sapphire ring laser (Co- 
herent MBR-110) operating at approximately 960 nm is 



frequency doubled in an external ring resonator (Coher- 
ent MBD-200) to produce approximately 90 mW at 480 
nm. Light pulses of 1 iis duration are produced using an 
acousto-optic modulator. After the acousto-optic modu- 
lator the beam is coupled into a single mode fiber. The 
output beam from the fiber is coUimated then focused 
and directed towards the trapped atoms. The highest 
Rydberg densities (w 2 x 10^ cm"-^) in this paper were 
obtained with optical powers of approximately 25 mW. 

Excitation of cold atoms to Rydberg states occurs as 
a two-color, two-photon process with the 480 nm light, 
and the nearly resonant, red detuned 780 nm light used 
for cooling and trapping [l3|. Unfortunately, there are 
several atmospheric water absorption lines in the vicinity 
of 960 nm which restrict performance of the Ti:sapphire 
laser system. To avoid one of these, we excite 46^5/2 
states instead of the 45^5/2 states used in Ref. [13 • 

Microwave horns are aimed at the Rydberg atoms 
through a fused silica window. To drive the 46^5/2 — 
47d5/2 (« 35.7 GHz), the 46^/2 - 47p3/2 (« 22.1 GHz), 
and the 46(^5/2 — 45/5/2 (~ 24.1 GHz) transitions we 
directly use the unamplified output of Agilent E8254A 
microwave synthesizers. For the 46^5/2 — 46p3/2 (« 
100.1 GHz) and the 34si/2 - 34pi/2 (« 104.1 GHz) tran- 
sitions both an active quadrupler (Spacek Labs P/N 
A100-4XW) or a passive tripler (Pacific MiUimeter P/N 
W3W0) were used. These multipliers were driven by an 
Agilent E8254A synthesizer. 

Two non-magnetic stainless steel electrode plates, sep- 
arated by 36 mm, are located to either side of the trapped 
atoms. These electrodes are used for SFI of Rydberg 
atoms [13. The resultant ions are drawn by the same 
fields through a 6 mm hole in one of the electrode plates 
towards a microchannel plate detector. These plates are 
also useful for controlling the electric field in the experi- 
mental region. By varying the relative voltage difference 
of the two plates, their average voltage and the voltage 
of the Rb dispenser source with respect to the grounded 
chamber, it is possible to compensate for stray electric 
fields in all three directions in the experimental region. 
Stray electric fields are measured using the Stark shifts 
of a one-photon microwave transition |l8l | . The residual 
electric field is estimated to be less than 0.1 V/cm. 

To obtain estimates of the Rydberg density it is nec- 
essary to know the number of Rydberg atoms and their 
spatial extent. The spatial profile of the 480 nm beam 
determines the extent of the Rydberg cloud in two direc- 
tions (the focus of this beam is much smaller than the 
trapped atom dimensions). Beam characterization can- 
not easily be done within the vacuum chamber. Instead, 
the rigidly attached output coupler and focusing assem- 
bly are moved to an optical table and scanning knife 
edge beam profile measurements are performed through- 
out the focal region (see for example Ref. 19] ). The 
beam profile is Gaussian. Since there are collision pro- 
cesses which depend on the Rydberg density, the loca- 
tion of the tightest focus in the trap can be determined 
by moving the focusing lens relative to the trap to max- 



imize these processes. Once this location is found, the 
excitation beam profile can be systematically varied by 
moving the focusing lens. The beam was set to have a 
full width at half maximum (FWHM) of 0.190±0.015mm 
at the location of the trap. 

The extent of the Rydberg cloud in the third dimen- 
sion is determined by the trapped atom profile. By scan- 
ning the direction of the 480 nm excitation laser - moving 
the beam over the trapped atom cloud - it is observed 
that there is an excellent correlation between the Ryd- 
berg signal and 780nm fluoresence imaging measurements 
of cloud extent (accounting for the finite width of the 
480 nm beam). Therefore 780 nm fluorescence imaging 
gives the spatial extent of the Rydberg cloud in the di- 
rection of the 480 nm beam propagation. This is typically 
0.5 ±0.1 mm (FWHM). 

The beam profile of the 480 nm Rydberg excitation 
laser and the spatial profile of the trapped atoms are 
combined to convert the observed field ionization signals 
into spatially averaged Rydberg densities. If the waist of 
the excitation laser is significantly smaller than the width 
of the trap, and this waist docs not change appreciably 
as the beam passes through the trap (conditions that 
are approximately satisfied in the present setup) the spa- 
tially averaged Rydberg density corresponding to a total 
number of Rydberg atoms Nnyd is 
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where FWHMj^ is the full width at half maximum of the 
irradiance profile of the 480 nm light, and FWHM|| cor- 
responds to the 5^3/2 spatial distribution in the direction 
of the 480 nm beam. 

Determination of N^yd requires knowledge of the ab- 
solute detection efficiency and gain of the microchannel 
plate detector (MCP). The absolute detection efficiency 
of an MCP depends on the kinetic energy at impact, 
the species involved, and surface conditions [23|- A rec- 
ommended efficiency for our specific conditions was not 
available. Based on information from the manufacturer 
and the literature p(]| | we estimate an initial detection 
efficiency of 50% - however this could be as high as 80% 
or as low as 20%. The gain of the MCP was calibrated 
by observing the average total charge due to single events 
at higher bias voltages than normal operating conditions. 
To determine gain at the operating bias voltage, a series 
of measurements were taken at constant charged particle 
signal level with variable bias voltages. 

The average densities quoted in this paper are obtained 
using Eq. ^ The uncertainty in these is dominated by 
NRyd, which is limited by the initial detection efficiency 
estimate. 
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FIG. 1: a) Timing for experiment. There are delays of 
roughly 100 ns between the excitation laser pulse, transfer 
pulse, probe pulse and field ionizing pulse. This sequence is 
repeated at 10 Hz and is synchronized with the 60 Hz AC 
power line, b) Energy levels of relevant states. The spectro- 
scopic data is from Ref. |2lll . The matrix elements magnitudes 
|<n^j|r|n'^'/>| are shown between appropriate states, where 
r is the radial coordinate of the Rydberg electron. These 
are expressed in atomic units, and were calculated using the 
techniques described by Zimmerman et al. |2^. 



III. OBSERVATION OF THE RESONANT 
ELECTRIC DIPOLE-DIPOLE INTERACTION 



The interactions between ®^Rb Rydberg atoms were 
studied using the linewidths of the 46(j .c^/9 — 47^5/2 two- 
photon "probe" transition. Li et al. [21| have shown 
that two-photon transitions between states with identi- 
cal gj factors, such as 46d5/2 and 47^5/2, show negligi- 
ble broadening in the inhomogeneous magnetic field of 
a MOT. These transitions also have less sensitivity to 
electric fields than one-photon transitions. 

The typical experimental sequence for a single shot is 
shown in Fig.^). Depending on the particular situation, 
a "transfer pulse" may be used shortly after optical ex- 
citation to drive a fraction of the initially excited 46^5/2 
atoms into other states. In the presence of the inhomoge- 



neous magnetic field of the MOT any coherent superposi- 
tions formed by this pulse will rapidly dephase across the 
sample. Interatomic interactions are somewhat inhomo- 
geneous due to the variation in spacings between atoms, 
and are also expected to cause appreciable dephasing. As 
a consequence the transfer pulse is viewed as redistribut- 
ing population rather than creating superpositions. 

The 46^5/2 ~ 47(i5/2 probe microwave frequency is 
scanned between shots. The power of this probe is ad- 
justed so that power broadening is negligible, and in the 
limit of low Rydberg densities, the linewidths of this 
transition are entirely dominated by transform broad- 
ening. To analyze additional broadening due to inter- 
atomic effects, the spectra are fit to the transform lim- 
ited lineshape for a square excitation pulse (sine (tt/T)) 
convolved with a Lorentzian of variable width Sv. The 
widths of the fitted Lorentzian are then analyzed as 
a function of spatially averaged total Rydberg density 
(which is normally adjusted using the 480 nm excitation 
laser power). 

As illustrated in Fig. [2 61^ converges to zero as the Ry- 
dberg density goes to zero. This indicates that the line 
broadening is an interatomic effect. At 1 x 10^ cm^-^ Ry- 
dberg density, transferring one-half of the 46^5/2 atoms 
to 47p3/2 broadens the lincwidth, increasing Sv from ap- 
proximately 10 kHz to 120 kHz (the uncertainty in Sv is 
typically ±7kHz). This can be attributed to the resonant 
electric dipole-dipole interaction between atoms in these 
two states. Consider a pair of atoms labelled A and B. 
In the absence of interactions between these atoms, the 
state |46(i5/2TO5/2>A|47p3/2m3/2>B is energy degenerate 
with |47p3/2m3/2>A|46d5/2TO5/2>B- The dipole-dipole 
interaction breaks this degeneracy, leading to an energy 
splitting on the order of 
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where r is the separation of Rydberg electron from the 
ion-core, and -Rab is the typical interatomic spacing |23 . 

The value of | < 46^5/2^147^3/2 > | and other rele- 
vant radial transition dipole moments are indicated in 
Fig. n^). For instance, the dipole coupling between the 
47(i5/2 and 47^^3/2 states is relatively weak compared to 
the coupling between the 46^5/2 and 47p3/2 states. Thus 
the broadening of the two-photon 46^5/2 — 47^5/2 probe 
transition beyond the transform limit is dominated by 
dipole-dipole energy splittings in the initial state. As in 
Ref. [13 , a simple estimate based on Eq. [3 is consistent 
with the observed broadening. The l/i?AB scaling of the 
splittings implies a linear density dependence of the line 
broadening, which is clearly demonstrated in Fig.[2Ji). 

The dipole-dipole coupling between 46^5/2 and 45/ 
states is of comparable magnitude to that between 46^5/2 
and 47p3/2 states (see Fig.QJ)). Consequently, transfer of 
a fraction of atoms to the 45/ state gives similar broad- 
ening (see Fig. [33). 



It is also straightforward to drive transitions to a 
nearby state which is not strongly dipole coupled to ei- 
ther the initial or final state of the two-photon probe 
transition. Based on the radial dipole matrix elements 
we would expect 4 orders of magnitude less broadening 
if 50% of the atoms were transferred to 46p3/2 instead 
of 47p3/2- As Fig. [2b) indicates, there is no detectable 
broadening of the probe transition in the presence of 
46p3/2 atoms. 

These results are similar to those presented in Ref. [13 
(only at one higher principal quantum number). There 
is a slight increase in linewidth of w 20% in a mixture 
of nd^/2 E^nd (n -|- 1)^3/2 atoms (where n = 46 here, and 
n = 45 in Ref. [13 )■ The energy splittings Avdd scale 
like n* due to the n^ scaling of the transition dipole mo- 
ments |l^, and this accounts for about a 10 % increase. 
We have also made changes in our Rydberg density es- 
timate procedure, and these may also account for the 
discrepency. The more significant difference is in the 
linewidths obtained when no transfer pulse is present (see 
Fig.[2J). There is roughly 3 times less broadening as com- 
pared to the data presented in Ref. 13]. 

The source of broadening in the case of no transfer 
pulse was proposed to be due to the collision of hot Ryd- 
berg atoms (w 300 K) with the cold Rydberg atoms |13l |. 
These types of collisions have also been discussed in the 
context of cold plasma formation [23 . The hot atoms are 
excited in the same way as the cold trapped atoms. They 
make an undetectable contribution to the observed field 
ionization because the waiting period between excitation 
and extraction allows them to move to where they can't 
be detected. 

Several experimental tests support the role of the hot 
Rydberg atoms in line broadening. To control the hot 
background Rb number density, the current through an 
alkali metal dispenser source is varied |2g| . When the hot 
Rb background density is varied in this manner, the 480 
nm excitation laser power should be adjusted to produce 
the same cold Rydberg density (compensating for the 
different trapped atom density). The time constant is 
relatively long for changes in background pressure with 
dispenser current (on the order of 10 minutes), so we can 
check independently of background Rb pressure whether 
thermal radiation or charged particles from the dispenser 
are responsible for broadening. At our sensitivity, no 
such effects were found. 

With constant cold Rydberg density, it is found that 
5v increases with hot Rb background density. This rules 
out the interactions between cold Rydberg atoms as be- 
ing responsible for the broadening. However, at higher 
densities than those studied here, the nearly resonant 
interactions 46^5/246^5/2 <-> 48^3/244/ between cold Ry- 
dberg atoms will become increasingly important [l2| . At 
a fixed hot Rb background density, Sv is proportional 
to 480 nm excitation laser power, which rules out the 
influence of unexcited hot atoms. 

Not only does the hot Rydberg atom density in- 
crease with background Rb pressure, but so does the 
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FIG. 2; Broadening of the 46^5/2 — 47^5/2 probe transition 
as a function of average Rydberg density with and without 
transfer pulses. The transfer pulses create 50 % mixtures, 
except in the case of the 46d5/2 -I- 45/ mixture, where only 
30 % of the total Rydberg atoms are transferred to the 45/ 
state. The straight lines are least squares fits. See the text 
for the definition of 5u. 



trapped atom density. This suggests a check to deter- 
mine whether interactions with less excited cold atoms, 
or with hot Rydberg atoms from the vapor, are the cause 
of the increased broadening. If the cooling and trapping 
laser frequency is detuned to obtain less cold atoms, and 
the 480 nm laser power is increased to obtain the same 
cold Rydberg density, the observed linewidths increase, 
because more hot atoms are optically excited to Rydberg 
states. If less excited cold atoms were responsible for the 
broadening we would expect a decrease in linewidth. 

In summary, the broadening 61/ observed with no trans- 
fer pulse is due to the collision of hot Rydberg atoms 
with the cold Rydberg atoms we are studying. In the 
present work Sv in the absence of a transfer pulse is 
smaller than in Ref. [13 for at least two reasons. We 
are now working at lower dispenser currents (correspond- 
ing to lower background pressures) and using a 480 nm 
excitation source with significantly narrower bandwidth. 
The reduced bandwidth is expected to alter the ratio of 
cold to hot Rydberg atoms. 

By working with a low hot Rb background pressure, 
and at 10^ cm^"^ total Rydberg density, 5v is on the or- 
der of its experimental uncertainty ±7 kHz. We have not 
yet established what the ultimate limit on the lowest 6v 
is, and thus the ultimate sensitivity of the probe tran- 
sition to Rydberg interactions. However, the currently 
demonstrated sensitivity is vastly superior to what can 
readily be achieved with optical probes. 



IV. MAGNETIC FIELD MEASUREMENTS 

In this section we discuss the use of a single-photon 
microwave transition between Rydberg states to measure 
magnetic fields. 

The MOT used to trap the ^^Rb atoms requires a mag- 
netic quadrupole field |l5| . In these experiments the field 
gradients in the three orthogonal directions are estimated 
to be 12, 12 and 24 G/cm, based on the coil geometry. 
The Rydberg cloud has a Gaussian profile with FWHM's 
of 0.190 ± 0.015 mm in two dimensions (dictated by the 
480 nm excitation laser) and 0.5 ±0.1 mm in the third 
dimension (dictated by the trap size) . The largest dimen- 
sion corresponds to the direction of the largest magnetic 
field gradient. This gives a root mean square (RMS) mag- 
netic field of approximately 0.5 G. Since Zeeman shifts 
are typically on the order of 1 MHz/G, and the dipole- 
dipole interactions that we observe are roughly 100 kHz, 
it is apparent that the quadrupole field is a significant in- 
homogeneity, and it is desirable to turn this off. The cur- 
rent in the coils producing this field is shut off in 1.7 ms, 
but significant eddy currents persist for much longer (see, 
for example, Ref. [23 )• Excitation to Rydberg states oc- 
curs after a delay of 25 ms, to allow the eddy currents to 
die down (see Fig. ^. As will be discussed, an upper 
bound can be placed on any remaining inhomogeneity. 

The detuning and counterpropagating beam balance of 
the MOT lasers are optimized to minimize the reduction 
in cold atom density during the 25 ms delay. Immediately 
after selective field ionization is complete, the current to 
the quadrupole coils is turned back on to recapture cold 
atoms for the next shot. 

Once the quadrupole field is turned off, a DC magnetic 
field may remain. Three pairs of orthogonal Helmholtz 
coils are used to null any preexisting field (due to the ion- 
pump, earth etc..) and one of these coil pairs is used to 
produce a deliberate homogeneous field. The current in 
the coils is constant throughout the MOT collection, Ry- 
dberg excitation, and microwave probe cycle. When set- 
ting a non-zero homogeneous field, the quadrupole zero 
of the MOT will be shifted. As long as the quadrupole 
zero remains significantly within the overlap of the cool- 
ing and trapping beams, the trap will still function. The 
use of fields of up to 0.6 G pointing in the direction with 
a field gradient of 12 G/cm, shifts the zero by 0.5 mm, 
which does not significantly influence MOT operation 
when using 780 nm beams with a FWHM of 4 mm. Loss 
of the cold atoms in the 25 ms waiting period limits the 
DC bias fields which can be applied using this technique. 

Although it is straightforward to calculate the fields 
due to the bias coils from their geometries, we do not a 
priori know the stray magnetic field within the appara- 
tus. To determine this, the one-photon microwave transi- 
tion 34si/2 — 34pi/2 of ^^Rb is used. We work at a lower n 
compared to where the dipole-dipole experiment is per- 
formed so that the hyperfine spacings are larger (these 
scale like 1/n^, see Ref. [13). To improve the sensitivity 
of the measurement it is desirable to use long microwave 



10q 

5-: 
0- 



Rydberg 
excitation 



E 1 
o 

O 

X 

r^ 

CD 

ro 

0. 
0. 
cq" 0. 

0. 



2-1 
8- 

4- 
0- 

6-1 
4- 
2- 
0- 



n-<- 

-30 



-25 



' I " " I " " I ' 

-20 -15 -10 
time (ms) 



FIG. 3: Timing for experiments done in the presence of a 
homogeneous DC magnetic field. The quadrupole field coil 
current is cut to zero in 1.7 ms. As the magnetic field gradi- 
ent dBx/dx decays to zero, the field in the center of the trap 
Bx gradually increases due to the presence of a non-zero ho- 
mogeneous bias field. At t = 0, excitation to Rydberg states 
and subsequent microwave probing takes place, as shown in 
Fig. The gradient and bias field time constant was mea- 
sured using Rydberg atom microwave spectroscopy (see text). 



pulse lengths (for narrow linewidths). Redistribution of 
Rydberg state population due to thermal radiation, lim- 
its the microwave pulse length. Pulses of 36 /iS are used. 

The Zeeman splittings of the relevant levels are shown 
in Fig. ^ The 780 nm light used for optical excitation is 
nearly resonant with only the 5si/2i^ = 3 to 5p3/2-F — 4 
transition. Since F cannot change by more than one in a 
single photon transition, the 34si/2F = 2 state is not ex- 
cited in zero magnetic field, and only two microwave spec- 
tral lines are observable: 34si/2,^ = 3 to 34:pi/2,F = 2 
and 34si/2, F = 3 to 34pi/2, F = 3 (see Fig. 0. 

As a magnetic field is applied, this spectrum splits into 
many lines. To zero the magnetic field in a specific di- 
rection the following procedure is followed. A current is 
applied to a set of coils to produce a relatively large mag- 
netic field (the precise value is unknown at this point). 
By scanning the microwave frequency, a spectrum with 
many lines can be observed (see Fig.O. Then we deter- 
mine the current which produces the same field magni- 
tude in the opposite direction, by trying to reproduce the 
spectral line positions. The average of these two currents 
corresponds to zero field projection along this axis. Us- 
ing this procedure all three orthogonal axes can be zeroed 
independently. Since the change in magnetic field with 
current can be estimated from coil geometry, the pre- 
cision to which the average currents can be determined 
can be directly related to the precision to which the fields 
can be compensated. We have zeroed the magnetic field 
along each axis to within ±6 mG using this technique. 
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FIG. 4: Calculated Zeeman structure of the 34si/2 and 
34pi/2 levels of *^Rb using the Breit-Rabi formulae |2a| (the 
mj labelling is approximate). Please note the different ver- 
tical scalings. The two hyperfine splittings are chosen to 
match the spectral lines in Fig. |3 (0.567 MHz for 34si/2, and 
0.147MHz for 34pi/2). The arrow labels the 34si/2(mF = -3) 
to 34pi/2('TiF = — 3) line, which is used for magnetic field cal- 
ibration (see text). 



To deliberately create magnetic fields in a specific di- 
rection, one of the coil currents was adjusted away from 
its "zero" value. The change in magnetic field with coil 
current can be estimated from the coil geometry. How- 
ever, the line positions in the 34si/2 — 34pi/2 spectra (see 
Fig. EJ can be used to directly measure the field. The 
shifting of energies is given by the Breit-Rabi formulae 
(see, for example Ref. |28|). In general, line positions are 
dependent on the magnetic field, and the hyperfine split- 
ting of both the 34si/2 and 34pi/2 states. The 34si/2 
hyperfine spacing is not known a priori. However, there 
is a particular line (marked with arrows in Fig. 0) , cor- 
responding to the transition from 34:Si/2i'mF = —3) to 
34pi/2(f7ii? = —3) whose shift in position with magnetic 
field does not depend on the hyperfine spacing. The tran- 
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FIG. 5: The Rb 34si/2 to 34pi/2 microwave spectra taken 
with different DC magnetic fields (top 7 traces) and with the 
MOT quadrupole field present (bottom trace). Spectra are 
offset vertically for clarity. The calculated line positions of 
the Amp = transitions using the energy levels from Fig. |1] 
are shown for the 0.65 G data (the microwave polarization 
vector is in same direction as the applied magnetic field). The 
arrow labels the 34si/2 (jtif = —3) to 34pi/2("T-F = —3) line, 
which is used for magnetic field calibration (see text). 



sition energy for this line is (2/3)fiBB, where fiB is the 
Bohr magneton and B is the magnetic field. As apparent 
in Fig. [SI this line is convenient as it is readily identifable 
in the spectra. By tabulating the position of this spec- 
tral line as a function of coil current and fitting a straight 
line, the magnetic field as a function of coil current can 
be determined, without knowledge of the 34si/2 hyper- 
fine splitting. This provides a precise determination of 
the magnetic field (±10 mG). The dipole-dipole interac- 
tion is then studied under these calibrated magnetic field 
conditions. 

It is desirable to know the inhomogeneity in the mag- 
netic field when the experiment is performed. The bot- 
tom trace in Fig.|Slshows inhomogeneous broadening due 
to the quadrupole field when the coil current is not shut 
off. By varying the time between between shutting the 
current off and excitation of the atoms to Rydberg states 
(see Fig. O we can observe the decay of this inhomogene- 
ity using the spectral widths. The decay is observed to 
be exponential with a time constant of 6 ± 1 ms. A gauss- 



meter on the outside of the vacuum chamber, placed as 
close as possible to the experimental region (within 15cm) 
shows a similar time constant for the change in magnetic 
field: 6.5 ± 0.5 ms. 

At the start of this section, the RMS magnetic field 
due to the quadrupole field was calculated to be 0.5 G. 
Based on the exponential decay of the gradient, 25 ms 
after the coil currents are shut-off we would expect this 
to decay to 8 mG. This inhomogeneity is insignificant 
in the dipolc-dipole interaction experiments discussed in 
the next section. 

The widths of the resolved 34si/2(?71f = —3) to 
34pi/2(f7iF = —3) spectral lines in Fig. |S1 put an upper 
bound on the magnetic field inhomogeneity. The lines 
are observed to have a linewidth of 30 kHz. The trans- 
form limited linewidth for a 36 ^s long pulse is approx- 
imately 25 kHz. If as a worst case we assume that the 
linewidth contributions combine in quadrature, the resid- 
ual broadening due to magnetic field inhomogeneity and 
other mechanisms is 17 kHz. This gives a rough upper 
bound of 17 mG on the magnetic field inhomogeneity, 
which is consistent with the estimate given in the previ- 
ous paragraph. 



V. SUPPRESSION OF THE RESONANT 

ELECTRIC DIPOLE-DIPOLE INTERACTION 

DUE TO A MAGNETIC FIELD 



Using the techniques described in the previous section, 
cold Rydberg atom interactions can be studied in the 
presence of homogeneous, calibrated magnetic fields. The 
timing of the experiment remains the same as in Fig.^i)- 
However, the current to the quadrupole coils is shut-off 
25 ms prior to photoexcitation (see Fig. (SJ. During this 
waiting period, the cold atom cloud dimensions expand. 
At each different applied magnetic field value, the coun- 
terpropagating beam powers are re-optimized to mini- 
mize expansion. When this is done, it is found that trap 
expansion is independent of DC magnetic field to within 
10% for fields less than 0.6 G, and linewidth broadening 
can be studied at a fixed Rydberg density as a function 
of magnetic field. 

The linewidth broadening due to interatomic effects, 
bv is obtained in exactly the same manner as described 
in Section imi In the limit of low densities bv approaches 
zero for all DG magnetic fields, indicating that it is due to 
interatomic interactions. At an average Rydberg density 
of 8.0 X 10^ cm"'^ we have found that conversion of 50 % 
of the 46(^5/2 atoms to 47p3/2 increases bv from 10 kHz 
to 110 kHz at approximately 0.04 G. As Fig. indicates, 
increasing the magnetic field reduces bv. The presence of 
a DG magnetic field suppresses the resonant dipole-dipole 
interaction between atoms. Essentially the magnetic field 
spoils some of the unperturbed atom energy degeneracies, 
weakening the resonant interaction. 

To discuss this in detail, we consider the dipole-dipole 
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FIG. 6: Broadening of the 46^5/2 — 47d5/2 probe transition in 
the presence of 50% 47p3/2 atoms as a function of magnetic 
field strength. Two diflterent total Rydberg number densi- 
ties are shown. The solid lines are calculations (see text for 
details) . 



interaction operator: 



Vdd^ 



flA ■ MB - 3(^A • n){flB ■ n) 



-"ab 



(3) 



where jIa and jIb are the electric dipole matrix element 
operators evaluated on each atom, n is the unit vector 
pointing between the atoms, and -Rab is the separation 
of the two atoms. 

The tensor product of the separated, non-interacting 
atom energy eigenstates can be used as a basis set for 
the evaluation of Vdd- The axis of quantization for the 
projection of the angular momentum of each atom is cho- 
sen to point in the same direction as the magnetic field. 
This is in contrast to the more conventional choice for 
interatomic interactions, where the quantization axis is 
chosen along n (see, for example, Ref. |23|). This choice 
complicates the evaluation of the matrix elements of Vdd, 
but simplifies the evaluation of the Zeeman part of the 
Hamiltonian. 

If we label each atom as A and B, in zero magnetic field 
with no dipole-dipole interaction, all 24 states of the form 
l>=|46d5/2TOj,Ai>A|47p3/2m-j,Bi>B are energy degener- 
ate with the 24 states of the form |2> = |47p3/2"T-j,A2>A 
46d5/2mj.B2>B. It is this degeneracy which is responsi- 
ble for the resonant electric dipole-dipole line broadening. 

Magnetic fields lift some of this degeneracy through 
the Zeeman effect. If \^pb> is a basis set vector of the 
form given above, then: 



<i>b\Vz\ipb>^ {gjAmjA + ffjBn2jB)Ms-B 



(4) 



where gj is given by the standard formula {gj — 4/3 for 
the ^P3/2 states, and gj = 6/5 for the ^^5/2 states [28i] ) . 
In this basis set, with the axis of quantization along the 
direction of the magnetic field, there are no off-diagonal 
elements of Vz- 
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TABLE I: Non-zero matrix elements <2|Vdd|l> where |1>= 
|46d5/2mj=l/2>A|47p3/2mj=l/2>B and |2>=|47p3/2mj,A>A 
|46d5/27Tij,B>B, and i?AB ~ 28.5 fim. The axis for angular 
momentum quantization is in the direction of B (see text). 
The dipole moments jIa and jIb in Eq.|21are evaluated using 
the techniques described in Ref. ^22]. The Zeeman energy 
shifts of the |2> states relative to the |1> state are: /^mag,2 — 
Mmag,i = [<2\Vz\2> - <l\Vz\l>]/B (calculated using Eq.H. 



|2> 


<2|l>dd|l> (kHz) 


/^mag,2 /^inag,l 


rUj^A 


mj,B 


rrij^tot 


n II B 


n± B 


(kHz/G) 


-1/2 


-1/2 


-1 





55 


-3547 


3/2 


-1/2 


1 


21 


-11 


187 


1/2 


1/2 


1 


-147 


74 





-1/2 


3/2 


1 


52 


-26 


-187 


3/2 


3/2 


2 





45 


3547 



Once a magnetic field is applied only pairs of states for 
which mj,Ai = 'mj.,B2 and ruj^Bi = "t-j,A2, are degenerate. 
As an example, consider the state |46(i5/2TOj = 1/2 >a 

47p3/2TOj = 1/2 >B- This is coupled by Vdd to several 
states, as shown in Table IH These states are all energy 
degenerate in the absence of a magnetic field and interac- 
tions. Once a magnetic field is applied only one coupled 
state has an identical Zeeman shift and remains degen- 
erate: |47p3/2"ij = l/2>A|46d5/2Wj — 1/2>B. Since the 
application of a magnetic field lowers the energy degen- 
eracy, we expect less line broadening since some of the 
dipole-dipole couplings are no longer resonant. 

A quantitative treatment of the resulting linewidth 
suppression will now be considered. The essence of this 
treatment is the numerical diagonalization of a Hamilto- 
nian using a basis set of two atom states - the 24-1-24 = 48 
states mentioned above. All of these 48 states are energy 
degenerate in the absence of a magnetic field and dipole- 
dipole interactions, so the simplified Hamiltonian is 



H = Vz + Vm- 



(5) 



This corresponds to a 48 by 48 matrix for fixed values of 
n, i?AB and B. The matrix elements of Vdd require the 
transition dipole moments /Ia and /Ib (see Eq. (S)), which 
can be evaluated using numerical integration of the Ryd- 
berg electron wave- functions ;22] . The matrix elements of 
Vz are obtained from Eq. 0] The resulting Hamiltonian 
matrix may be numerically diagonalized, yielding a set 
of energy eigenvalues Ei and corresponding eigenvectors 

As has been discussed, the two-photon probe transi- 
tion is not sensitive to the Zeeman effect |2l[. Thus for 
a correspondence with the spectral transition energies, 
the Zeeman energy contribution is subtracted from each 
eigenvalue: 

E,, ^ E,- <^,\Vz\i^,> . (6) 

This heuristic expression is based on the observation that 



the probe linewidth is independent of magnetic field in 
the absence of interatomic effects. An exact treatment 
would take into account relative linestrengths in the two- 
photon probe transition. 

The corrected energies Ed may be computed as a func- 
tion of magnetic field magnitude for a fixed i?AB and 
relative orientation of n and B. The distribution of the 
corrected energies dictates the linewidth broadening. If 
there was no dipole-dipole interaction these would all be 
zero. The spread in the corrected energies decreases with 
magnetic field - the magnetic field suppresses the dipole- 
dipole interaction. To quantify this with a single number, 
we have computed the root mean square (RMS) corrected 
energy <E^^>^^^ as a function of magnetic field for n ± i? 
and ri || i? at a fixed value of Rab (see Fig. [TJ. 

The quantitative difference in the suppression for n _L 
B and n || -B may be understood from the difference in 
the selection rules for Vdd in these two cases. When n 
is parallel to the quantization axis (-B), the dipole-dipole 
interaction Vdd does not couple states of different fTi^.tot , 
where m^tot — 'mj,A + 'mj,B- This is a consequence of the 
invariance of Vdd to rotations about n. However the same 
selection rule does not exist when n is perpendicular to 
the quantization axis (see, for example, Table P). Due 
to the similarity in gj factors for the ^P3/2 and ^^5/2 
states (4/3 vs. 6/5) the Zeeman shifts for states of the 
same mj^tot are quite similar when compared to states of 
different TOj^tot (see the last column of Tabled). Therefore 
the degeneracy between states coupled by Vdd with n\\B 
is spoiled at much higher fields than for n ± B. If Fig. [7| 
is extended to higher fields, then <E^^>^/^ for n \\ B 
eventually drops to 67 kHz. 

In the sample both n and i?AB are distributed ran- 
domly. The nearest neighbor probability distribution is 
given by |30l| : 



dPjRAB) 
di?AB 



= 47ri?AB"-int exp(-47ri?^Bnint/3). (7) 



The number density n-mt in this expression is half of 
the total Rydberg density ~ as half of the atoms reso- 
nantly interact with the other half in an even mixture 
of 46^5/2 and 47^3/2 atoms. Using this expression and 
random orientations of n, we have calculated the prob- 
ability distribution of the corrected eigenvalues. These 
distributions were convolved with a sine (tt/T) lineshape 
(T = 6 fis) and fitted in the same manner as the exper- 
imental data (using the sinc^(7r/T) lineshape convolved 
with a Lorentzian of adjustable width Sv). 

To make a comparison with the experimental observa- 
tions, the calculations must be done at a specific number 
density. As discussed, our experimental density determi- 
nation is uncertain by almost a factor of 2. Therefore 
the density used in the calculations is adjusted for the 
best fit to the data. This is the only adjustable parame- 
ter in the calculations, and we use the same adjustment 
factor for both densities in Fig. El Specifically, at the 
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FIG. 7: The calculated RMS of the corrected energy eigen- 
values Eci due to the electric dipole-dipole interaction be- 
tween 46d5/2 and 47p3/2 atoms, as a function of magnetic 
field strength (see text). The interatomic separation is -Rab ~ 
28.5 fim, the most probable resonantly interacting neighbor 
distance at a total Rydberg density of 1.4 x 10^ cm~^. Calcu- 
lations are shown for two different relative orientations of the 
magnetic field B, and the unit vector pointing between the 
two atoms n. 



higher density the calculations are performed at a con- 
stant Rydberg density of T4 x 10^ cm~'^, whereas the 
spatially averaged experimental density is estimated to 
be (0.8 ± 0.5) X 10^ cm"-^. The discrepency may be an 
indication of the importance of multi-body interactions, 
which are expected to enhance the interactions over a 
simple binary model. 

The calculations shown in Fig.Elpredict a reduction in 
linewidths with magnetic field comparable to the obser- 
vations. However the agreement is not perfect. In partic- 
ular, the experimental results appear to extrapolate to a 
lower linewidth at zero magnetic field. There are several 
possible reasons for this. Hyperfine structure has been 
neglected. This is expected to be especially important at 
low magnetic fields where the Zeeman shifts are smaller 
than the hyperfine splitting. In addition, stray electric 
fields break degeneracies between states of different \mj\, 
and this is also expected to suppress linewidth. 

The distribution of energy eigenvalues of the two 
atom system has been used to obtain linewidths - but 
linestrengths for the different transitions should be ac- 
counted for. In addition, not all of the eigenstates will 
be equally populated by the excitation scheme. 

Perhaps the most significant problem with the theoret- 
ical model is the assumption of binary interactions: only 
pairs of atoms interact. Experimental observations of res- 



onant energy transfer between cold Rydberg atoms have 
demonstrated that this picture is not complete 0, UM- 
Taking these into account would require a much more 
sophisticated approach |j, IsJ . 

VI. CONCLUDING REMARKS 



The magnetic field induced partial suppression of the 
resonant electric dipole-dipole interaction between cold 
Rydberg atoms has been observed, and a quantitative 
model of this phenomena has been presented. This model 
could be tested in more detail. A "one-dimensional sam- 
ple" has recently been successfully employed to determine 
the influence of the relative orientation of an electric field 
with fi on Rydberg dipole-dipole interactions [sj • By fo- 
cusing the Rydberg state excitation laser more tightly, so 
that the average interatomic separation is larger than the 
beam waist, the sample becomes one dimensional _32|. 
This gives a preferred direction to n in the lab-frame, and 
thus makes it possible to study interactions with B Jl n 
or B 11 n. As Fig.[Z|indicates, there should be a dramatic 
difference between these two cases, which would provide 
a useful check. 

Homogeneous DC magnetic fields cannot completely 
remove the energy degeneracies responsible for the reso- 
nant dipole-dipole interaction. For example, in the ab- 
sence of interactions |46d5/2'7i5/2 >a |47p3/2"^3/2 >b is 
energy degenerate with |47p3/2m3/2 >a \^&d^/2'm^/2 >b 
independent of magnetic field. For complete suppression 
of the dipole-dipole interaction it is interesting to con- 
sider the use of magnetic field gradients. By making the 
energy levels of the atoms depend on position, the reso- 
nance conditions can be completely spoiled. Strong mag- 
netic field gradients have recently been used in electron 
spin resonance (ESR) to suppress spin diffusion due to 
the magnetic dipolar flip-flop interaction 33]. An anal- 
ogous experiment could be performed in the case of res- 
onant electric dipole-dipole interactions. The required 
field gradient is modest. For example, in the present 
work we typically have 100 kHz interaction strengths at 
interatomic separations of 30 ^m. With Zeeman shifts 
typically being on the order of 1 MHz/G, field gradients 
of roughly 50 G/cm should give significant supression. 
These gradients are readily achievable - in magnetic mi- 
crotraps for example i34i] . Strong magnetic field gradients 
could be useful for suppression of resonant dipole-dipole 
interactions in situations where they are undesirable. 
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